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Abstract. We investigate dark matter halo properties as a function of a time- varying 
dark energy equation of state. The dynamics of the collapse of the halo is governed 
by the form of the quintessence potential, the time evolution of its equation of state, 
the initial conditions of the field and its homogeneity nature in the highly non-linear 
regime. These have a direct impact on the turnaround, virialisation and collapse 
times, altering in consequence the non-linear density contrast and virial radius. We 
compute halo concentrations using the Eke, Navarro & Steinmetz algorithm, examining 
two extreme scenarios: first, we assume that the quintessence field does not exhibit 
fluctuations on cluster scales and below - homogeneous fluid; second, we assume that 
the field inside the overdensity collapses along with the dark matter - inhomogeneous 
fluid. The Eke, Navarro & Steinmetz prescription reveals, in general, higher halo 
concentrations in inhomogeneous dark energy models than in their homogeneous 
equivalents. Halo concentrations appear to be controlled by both changes in formation 
epochs of the halo cores as well as by differing virialisation overdensities. We derive 
physical halo properties in all models and discuss their observational implications. We 
examine two possible methods for comparing observations with theoretical predictions. 
The first method works on galaxy cluster scales and consists of fitting the observed 
X-ray cluster gas density distributions to those predicted for an NEW profile. The 
second method works on galaxy scales and involves the observational measurement of 
the so-called central density parameter. 
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1. Introduction 

Observations of the Cosmic Microwave Background temperature anisotropy reveal that 
a mysterious constituent with negative pressure, so-called dark energy, accounts for 70 
percent of today's mass-energy budget and is causing the expansion of the universe to 
accelerate [1]. These observations are in remarkable concord with the observations of 
distant supernovae [2]. 

The present-day challenge in cosmology is to discover the physical nature of 
dark energy. Candidates for this unknown entity include the cosmological constant, 
a dynamic cosmic field, commonly designated quintessence, which varies across space 
and changes with time, and a possible modification of General Relativity at large scales 
[3-5]. 

Although the physical nature of dark energy is undisclosed, one can explore its 
effects on cosmic structure formation, in particular one can study its implications for 
the number density of dark matter halos and their density profiles. In this respect, 
significant progress has been made by several authors who performed numerical N- 
body simulations in dark energy models [6-10]. These investigations are imperative for 
cosmological studies that rely on these ingredients to measure dark energy. Examples of 
this studies include semi-analytical studies of strong lensing statistics [11-13] and weak 
lensing number counts [11,14]. 

In this paper, we investigate how halo properties change in cosmological models 
with dynamical dark energy. This work extends upon previous studies in that we 
examine halo properties as a function of a time-varying dark energy equation of state, 
covering four types of potentials, and its homogeneity nature in the highly non-linear 
regime. We utilise the predictions of the spherical collapse model, such as the virial 
overdensity, obtained by [15,16], and calculate halo concentrations using the semi- 
analytical algorithm of Eke, Navarro & Stcinmctz (ENS) [17]. We then derive physical 
halo properties in all models and discuss their observational implications. 

The behaviour of linear perturbations in a scalar field and its effect on structure 
formation has been investigated by a number of authors. However, the behaviour of 
quintessence during the non-linear gravitational collapse is not well understood and 
is currently under investigation (see e.g. [15,16,23-25] for recent work). Usually, it is 
assumed that the quintessence field does not exhibit density fiuctuations on cluster scales 
and below. The reason for this assumption is that, according to linear perturbation 
theory, the mass of the field is very small (the associated wavelength of the particle is 
of the order of the Hubble radius) and, hence, it does not feel matter overdensities of 
the size of tenth of a Mpc or smaller [26]. 

The assumption of neglecting the effects of matter perturbations on the evolution 
of dark energy at small scales is indeed a good approximation when perturbations in 
the metric are very small. However, care must be taken when extrapolating the small- 
scale hnear-regime results to the highly non-hnear regime. Then, locally the fiat FRW 
metric is no longer a good approximation to describe the geometry of overdense regions. 
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Highly non-linear matter perturbations could, in principle, modify the evohition of 
perturbations in dark energy considerably, and these could, in turn, backreact and 
affect the evolution of matter overdensities. Moreover, it is natural to think that once 
a dark matter overdensity decouples from the background expansion and collapses, the 
field inside the cluster feels the gravitational potential inside the overdensity and its 
evolution will be different from the background evolution. This is a general feature 
of many cosmological scalar fields whose properties depend on the local density of the 
region they "live in" [43] . 

[35] suggested that the quintessence field could have an important impact in the 
highly non — linear regime. [24, 31] noted that the quintessence field could indeed be 
important on galactic scales. It was put forward by [30] that it could in fact be 
responsible for the observed fiat rotation curves in galaxies. Other authors [37] discussed 
more exotic models, based on tachyon fields, and argued that the equation of state is 
scale-dependent . 

If it turns out that the effects of dark matter density perturbations and metric 
influence perturbations of quintessence on small scales, this could significantly change 
our understanding of structure formation on galactic and cluster scales. [15, 16, 39] 
have shown that properties of halos, such as the density contrast and the virial radius, 
depend critically on the form of the potential, the initial conditions of the field, the 
time evolution of its equation of state and on the behaviour of quintessence in highly 
non-linear regions. In reality, the dependence on the inhomogeneity of dark energy is 
only important for some dark energy candidates. If the dark energy equation of state is 
constant, the differences between the homogeneous and inhomogeneous cases are small, 
as long as the equation of state w does not differ largely from w — —1 [15, 16]. Thus, 
for constant equation of state, the fitting formulae for the cold dark matter (CDM) 
density contrast presented in the hterature [26,27], do not change drastically, even if 
inhomogeneities in the dark energy component are taken into account. 

The paper is organised as follows. In section 2 we describe briefly the spherical 
collapse model and its dependence on the homogeneity nature of dark energy. In section 
3 we calculate halo concentrations in both homogeneous and inhomogeneous dark energy 
models. We diagnose observational methods to measure physical halo properties in 
section 4. Section 5 discusses the results and draws the conclusions. 

2. The spherical collapse model and the homogeneity nature of dark energy 

We consider a fiat, homogeneous and isotropic background universe with scale factor 
a{t). Since we are interested in the matter dominated epoch, when structure formation 
starts, we fill the universe with cold dark matter of density pm cx a^^ and a dark energy 
fluid with energy density p^. The equations that describe the background universe are 
(we set h — c = 1 throughout the paper): 



(1) 
(2) 
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Figure 1. Dark energy equation of state in the background universe as a function of 
rcdshift, z. Thick soUd hnc - [33], dashed hnc - [45], dashed-dotted Une - [29], dotted 
hne - [36], thin solid hne - cosmological constant. 



where H = a/a is the Hubble rate. When = —1 dark energy is the vacuum 
energy density. If dark energy is a scalar field (quintessence), p<^ = |</)^ + 1^(0) 
and Pff, = — V{(j)), where V{4>) is the scalar field potential. In this case, it is useful 
to rewrite equation (2) as 

+ 3H<j) + V' = (3) 

where the prime represents a derivative with respect to 0. 

In this paper, we consider four examples of quintessential potentials: 

• The double exponential potential [33]: 

V{(P) = M (exp(/?0) + exp(70)) . (4) 

• The exponential potential with inverse power [45] : 

y(0) = M(exp(7/0)-l) (5) 

• The Albrecht-Skordis potential [29] : 

V{^) = M (A + (0 - Bf) exp(-70) (6) 

• The supergravity-motivated potential [36]: 

V^(0) = Mexp(02)/0T (7) 

We choose the parameters in the potentials and the initial conditions in the background 
such that we obtain the following present-day cosmological parameters: Q^n = 0.3, 
= 0.7, Ho = lOOhkms-1 Mpc^^ with /i = 0.7, and -I < W4, < -0.8. The background 
time evolution of the equation of state for each model is shown in figure 1. 

The evolution of a spherical overdense patch of scale radius R{t) is given by the 
Raychaudhuri equation: 

3R = -AttGR (p^, + + 3w^J) (8) 
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Note that it would be wrong to use the Priedmann equation for a closed universe with 
a constant curvature k, since the former can vary in time [26,27]. In the halo, the 
evolution of and is given by 

R 

= -3-^(1 + w^c)P0c + r (9) 

and oc R^'^ due to mass conservation. Once again, in the case of a scalar field, 
equation (9) can be written as [15, 16] 

R 

'r 



^c + sl^ + KW^^ (10) 



c 



where 0^ is the field inside the overdensity (the lowerscript distinguishes it from the 
background value) and Vc = V{(f)c) is its potential. The quantity F describes the energy 
loss of dark energy inside the dark matter overdensity, as this component does not 
necessarily follow the collapse of dark matter, and energy can formally flow out of the 
system. As such, F encodes in how far dark matter non-linearities act on the dark 
energy component. 

We will make two assumptions for F. In the first case, the quintessence field is 
assumed to be smooth throughout space. This corresponds to 



r = -3 --^ 0^, (11) 




0c(^i) = 4>{'ti) ^-iid (pciti) = <i^{'ti), which implies (pdt) = 0(^) at all times, and is the 
standard assumption made in the literature. In the second case, we assume that the 
field follows the dark matter collapse from the very beginning. That is, we assume 
F = 0. 

Clearly, the values of F chosen are not reahstic, but they mark off two extreme 
scenarios. At very late times, during the collapse of the dark matter, especially when 
the density contrast in dark matter is very large (^m 3> 1), the field should no longer 
feel the background metric, i.e. it decouples from the background expansion. In this 
regime, the evolution of dark energy could be different and infiuence the details of the 
collapse. Only when F is quantified and the boundary conditions between the outer and 
inner metrics are properly understood, will the spherical collapse model be able make 
to solid predictions. An estimate of F can only be obtained from a general relativistic 
treatment. This could be achieved in two ways. One way would be to develop a Swiss 
cheese model. Another way would be to perform N-body simulations that take into 
account the effects of dark energy perturbations in smaller scales. This is however out 
of the scope of this article. 

Throughout the paper, we refer to a homogeneous dark energy model as a 
model where dark energy does not exhibit fluctuations on cluster scales and below 
(F given by equation 11). This designation may be seen as somewhat abusive 
in the sense that these models are inhomogeneous on scales larger than cluster 
scales. In truth, only the cosmological constant model is homogeneous across all 
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space. Wc refer to a inhomogeneous dark energy model as a model where dark 
energy does exhibit fluctuations on cluster scales and below (T = 0). The terms 
homogeneous/inhomogeneous are therefore associated with cluster scales and below, 
here taken as the scales of interest. 

We evolve the spherical overdensity from high redshift until its virialisation occurs 
(see [15,16] for details). Figure 2 plots the CDM density contrast, Ayir = as a 
function of redshift. It is clear from the figure that an homogeneous dark energy fiuid is 
very similar to a cosmological constant. Major differences occur in the inhomogeneous 
case. These are discussed in the following section. 



3. Halo properties 

The dark matter halo is modelled with the NFW profile [44]: 

where is the characteristic length and Ps is the characteristic density scale, Ps — ScPcrit, 
where pcrit = SH'^/SnG is the critical density for closure and 

<^c=^r^ , r (13) 

3 ln{l + Cvir) - Cvir/(1 + Cyir) 

is the characteristic density contrast. The halo concentration is defined as Cyir = t^u/ts 
where Tyir is the radius of a sphere containing a mean density Ayir times the critical 
density: 

47r 

Mvir = yr^i^AvirPcrit. (14) 

At small radii (r <^ r^) p(r) ~ r^^, and at large radii (r ^ r^) p(r) r^^ . The 
NFW profile appears to be a good fit to numerically simulated halos over a wide range 
of masses in various cosmological scenarios. Although the exact inner slope is under 
debate, there is general consensus that the density profile steepens at large radii. 



3.1. Halo concentrations 

Numerical simulations repeatedly show that the later a halo forms the less concentrated 
it is. This appears as a refiection of the smaller cosmic density at later cosmic epochs. 
In hierarchical structure formation models less massive halos form earlier than more 
massive ones, hence the former are more concentrated than the latter. 

We use the prescription of [17] to calculate dark matter halo concentrations. This 
prescription consists of a simple analytic algorithm that was found to reproduce well the 
mass and redshift dependence of concentrations obtained from high-resolution N-body 
simulations performed in various cosmological models, in particular in a fiat model with 
a cosmological constant. 

One should, however, point out that the ENS model is by no means the only 
viable prescription. In fact, it is just one of a number of theoretical models for how 
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Figure 2. left four inhomogeneous dark energy models (thick solid line - [45], dashed 
line - [36] , dot-dashed line - [29] , dotted line - [33] plus a cosmological constant model 
(thin solid line); right six homogeneous dark energy models (thick solid line - [45], thick 
dashed hne - [36], dot-dashed line - [29], dotted line - [33], thin dashed line - constant 
equation of state w = —0.6, thin solid line - cosmological constant); upper panel dark 
matter halo concentrations; upper middle panel non-linear overdensity at virialisation; 
lower middle panel non-linear overdensity at virialisation normalised at Zq', lower panel 
characteristic formation redshift as defined by [17]. 
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the concentration mass relation arises. Examples of other models are the original NFW 
model, the Bullock ct al. model [18], or the recent extension of this by Maccio et al. 
[19]. Some of these studies produced results in agreement with ENS (e.g. [21]), others 
find a mass and redshift dependence that contradicts the ENS model (e.g. [20], [22]). 
At present the matter has not been settled, and it is not possible to identify one model 
that works better than all of the others. In this work we will restric ourselves to the ENS 
model. A Further investigation using and comparing other prescriptions will appear in 
a near future work [50] . 

The ENS algorithm relates halo properties with the physics of halo formation as 
follows: 



which results from setting the core density to the spherical top-hat density at the 
characteristic formation epoch. Zo is the redshift at which the halo is identified which, 
is assumed to be the redshift at which the halo as a whole collapses and virialises. The 
characteristic formation epoch (also designated rapid-collapse epoch), Zc, is associated 
with the collapse time of the halo subunits. This occurs at the time when the rapid 
mass accretion rate drops below a fixed value. Thereafter, the halo evolves and its virial 
radius and halo mass grow through minor mergers and diffuse mass accretion, while the 
characteristic length remains essentially equal [8,47]. This process continues until 
the halo as a whole finally collapses and virialises. The rapid-collapse epoch of a halo 
of mass M depends on the linear growth factor, D{z), the amplitude and shape of the 
matter power spectrum, and on a single free parameter C„ whose value can be found 
which yields concentration values that match the results from N-body simulations. More 
explicitly, [17] define the collapse redshift as: 



where D{z) is the linear growth factor and aes is an effective amplitude of the power 
spectrum that [17] utilise to modulate a{M) in order to model WDM models. This 
effective amplitude is computed at Ms, the mass enclosed within the radius at which 
the NFW circular velocity reaches its maximum, r^ax — 2.17rs. 

[17] found that C^- = 28 gives good agreement with the N-body simulations 
effectuated in a fiat cosmological constant model. [7] have recently performed high- 
resolution numerical simulations in dark energy models with constant and time-varying 
equation of state (including Ratra-Pccblcs and SUGRA models) and have shown that 
the mass dependence and redshift evolution of concentrations is in conformity with the 
[17] prescription with C^. = 28. Moreover, these authors have demonstrated for the first 
time with numerical simulations that concentrations are indeed larger in dark energy 
models as a consequence of halos forming earlier in these models when compared to the 
cosmological constant model. This confirms the predictions of [34] and [12] which were 
based on the analytic algorithm of [17] 




(15) 



D{zc)(J^fi{Ms) 



1 



(16) 
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In this work we assume that the [17] prescription remains vahd for models 
with inhomogeneous dark energy, once all the cosmological functions inherent to the 
algorithm are specified accordingly. This may in fact be a strong assumption. In 
particular, the value assumed for the fitting parameter C„ = 28 might not be the 
correct one. However, this is the only reference value one has at present. Up until 
now, there are no numerical simulations which investigate the effects of dark energy on 
the internal dynamics of the halo formation. Dark energy has only been incorporated 
in the background evolution. Both its pressure and energy density have not even 
been considered to contribute to the local gravitational potential inside overdensities. 
Although in a very simplified way, the first studies which took into consideration such 
contributions were the N-body simulations performed by [41] in the case of coupled 
quintessence. 

For the shape of the linear matter power spectrum we adopt the fitting formula 
of [40] for QCDM models and the [32] fitting formula for the ACDM model using the 
modification of [46] in order to account for the baryons. We fix the primordial power 
spectrum index to n = 1 and the baryon density to fl^h'^ — 0.02. The normalisation 
of the power spectrum is set by as, the rms linear density fiuctuation in spheres of 
8h~^ Mpc, which we choose to be erg = 0.9. 

3.1.1. Inhomogeneous versus homogeneous models Figure 2 shows halo concentrations 
for all the models studied in this paper (on the left - inhomogeneous dark energy 
models plus a cosmological constant model for comparison, on the right - homogeneous 
dark energy models, including also the cosmological constant model for comparison). 
Also shown are the rapid-collapse redshifts, Zc (as defined by [17]) and the non-hnear 
overdensity at virialisation, Ayir, for the respective models. 

The upper panel indicates that the [17] prescription predicts, in general, larger 
halo concentrations in inhomogeneous dark energy models than in homogeneous ones. 
This is however model dependent. For instance, the [36] inhomogeneous model is the 
model with the highest concentrations which are nearly a factor of two higher than 
the cosmological constant model or its homogeneous counterpart. In contrast, the 
inhomogeneous [29] model is almost indistinguishable from the cosmological constant 
model and its homogeneous counterpart. The inhomogeneous [33] model and [45] model 
lie in between the cosmological constant model and the [36] inhomogeneous model. The 
former reveal concentrations higher than their homogeneous counterparts. Furthermore, 
all homogeneous models present concentrations of the same order of the cosmological 
constant model. 

This is interpreted as a combined effect arising from the ratios that enter as factors 
in formula (15). The matter density ratio is essentially equal in all the models. Although 
there is some difference on the rapid-collapse redshifts (lower panel), the overall change 
in concentrations for inhomogeneous models is primarily owing to the contrast in non- 
linear overdensity at viriahsation (middle panels). The latter depends on the clustering 
properties of the quintessence field. While at high viriahsation redshifts, all models 
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predict Avir ~ 178, significant deviations may occur at low virialisation redshifts due 
to the fact that dark energy starts to dominate in the background universe. Indeed, 
for inhomogeneous dark energy (F = 0), Avir can differ by a factor of four or more at 
low virialisation redshifts. This is nevertheless model dependent. In the [29] model, the 
field behaves like a cosmological constant in the background (see Fig 1). Hence, one 
would expect small differences between that and the ACDM model. This is indeed the 
case, because 0c ~ which, in turn, implies F 0, as it can be seen from equation (11). 
Thus, in this model, the fluctuations in the quintessence field remain small and the field 
is almost homogeneous. In fact, the non-hnear overdensity for the Albrecht & Skordis 
model oscillates around the cosmological constant model which explains the virtually 
coinciding concentrations in these two models. One should point out that what really 
characterises the differing ratios Avir/Avir(-2o), and consequent varying concentrations in 
the various models, is the value of the non-linear overdensity at Zo (here, Zo = 0.05) since 
at high redshifts, where the rapid-collapse of the halo happens, its value is practically 
equal in all models |. 

In effect, one can see that A^ir at low redshifts can differ by a factor of two between 
the cosmological constant model and an inhomogeneous model, but differs only about 
ten per cent between the cosmological constant model and an homogeneous dark energy 
model. As a result, homogeneous dark energy models exhibit less noticeable differences 
in concentrations (unless one considers high values of w) and these are mostly inherent 
to the differing rapid-collapse redshifts which arc highest for the model with constant 
equation of state w = —0.6. Indeed, for homogeneous dark energy models concentrations 
are larger as w increases owing to structures forming earlier. This is in agreement with 
the previous findings from high-resolution numerical simulations [6-8] which confirm 
that halos keep a memory in their central regions of the mean density of the universe 
at the characteristic formation epoch. Why do all the models, both homogeneous and 
inhomogeneous models, exhibit analogous values of Zc? It is because Zc is calculated 
from linear theory, depending solely on the matter power spectrum and on the linear 
growth factor. Hence, this depends on the background cosmology only and not on the 
non-linear physics dynamics of the halo that distinctively categorises the inhomogeneous 
dark energy models and from which one can infer the linear overdensity at collapse. The 
omission of the latter in the [17] prescription may be invalid, and thus needs to be tested. 

To summarise, the trends presented in Fig. 2 for halo concentrations in the two 
classes of models studied appear to be controlled both by differing formation histories 
as well as by varying virialisation overdensities. While for inhomogeneous models, the 
virialisation overdensities constitute the dominant factor, for homogeneous models it is 
the differing formation histories that account for the changes in concentrations. 
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Figure 3. left ps vs. rs relation in four inhomogeneous dark energy models (thick solid line - 
[45] , dashed line - [36] , dot-dashed line - [29] , dotted line - [33] plus a cosmological constant 
model (thin solid line); right the same relation in six homogeneous dark energy models (thick 
solid line - [45] , thick dashed line - [36] , dot-dashed line - [29] , dotted line - [33] , thin dashed 
line - constant equation of state w = —0.6, thin solid line - cosmological constant). 




Figure 4. left A^,/2 vs. V^max in four inhomogeneous dark energy models (thick solid line - 
[45], dashed line - [36], dot-dashed line - [29], dotted line - [33] plus a cosmological constant 
model (thin solid line); right the same relation in six homogeneous dark energy models (thick 
solid line - [45], thick dashed line - [36], dot-dashed line - [29], dotted line - [33], thin dashed 
line - constant equation of state w = —0.6, thin solid line - cosmological constant). 



4. Observational Implications 

The effect of dark energy on halo properties is expected to be imprinted on observable 
quantities. Unfortunately, the predicted Cyir — M^i,. relation is not itself directly 
observable. 

Note that [38] have determined the Cyir — Myij. relation using the dynamical 
properties of dark halos derived from the distribution of galaxies in massive systems 
(groups or clusters) and the rotation curves in less massive systems (dwarf, low surface 

I it approaches the Einstein-de-Sitter model Ayir ~ 178. 
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brightness and spiral galaxies) and have found that the observational data are grossly 
consistent with the theoretical predictions of a flat cosmological constant model on a 
mass range 10^° — 10^^ Mq. However, measurements of halo concentrations through the 
aforementioned observational method, as well as via weak and/or strong lensing cluster 
measurements, cannot be used to discriminate cosmological models since observational 
estimations of their values rely on Avir and Pcrit which are cosmology-dependent but 
unobservable. 

In the following, we examine two possible methods for comparing observations with 
theoretical predictions. The first method works on galaxy cluster scales and consists 
of fitting the observed X-ray cluster gas density distributions to those predicted for an 
NFW profile. [42] have analytically shown that the density distribution of an isothermal 
gas cloud with temperature Tx, in hydrostatic equilibrium, within an NFW dark matter 
halo is: 

Pgasir) = P,a.(0)e-"(1 + r/rj^/^''/^^) (17) 

where a — ATrG/impPsrg/kTx and /i and rrip designate the mean molecular weight 
and the proton mass, respectively. The best-fit parameters are thus a and r^, or, 
alternatively, if one utilises the X-ray temperature, ps and r^. [48] have applied this 
technique to an ensemble of 63 X-ray luminous clusters and have determined (assuming 
a fiat cosmological constant model) a best-fit Cvir — M^i^. relation. Here, we look directly 
at the quantities ps and which are not dependent on cosmology-dependent but 
unobservable quantities such as iiVir, although comparison with observations rests on 
assumptions such as the hydrostatic equilibrium and an isothermal temperature profile 
of the intracluster gas. 

The second method works on galaxy scales and involves the observational 
measurement of the so-called central density parameter proposed by [28]. This 
parameter is defined as the mean dark matter overdensity within the radius, rv/2, where 
the galaxy rotation curve is one half of its maximum, Vmax- 

The link between theory and observations can be made through the Ay/2~Knax relation, 
where Knax works as a measure of the absolute size of the halo. As pointed out by [49], 
Av/2 has the advantage of being defined without reference to any particular density or 
velocity profile. Theoretically, the maximum velocity of a NFW halo is given by [49]: 

V^^ ~ 0.2161^4 — (19) 

Z77;(Cvir) Cyir/ (l ~l~ Cyir) 

where Kir is the virial velocity and V^^^ = GMvir/-Rvir- In addition, rv/2 — 0.13rs 
for a NFW profile. Incorporating equation 19 into equation 18 yields the expected 
Ay/2 in theory. This can be computed for any cosmological model given Cyir and Ayir. 
Comparison with observational data relies, however, on the crucial assumption that 
baryons have not substantially modified the dark matter profile. This assumption is 
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weakened, but not unambiguously removed, if observations are focused on low-surface 
brightness and dwarf galaxies which are believed to be dominated by dark matter. 

We now look at the theoretical predictions from both methods. Fig. 3 plots the 
theoretical relation ps — rg expected in the various cosmological models for halos at 
redshift Zo = 0.05. Fig. 4 presents the theoretical relation Ay/2 — Knax expected 
in the same cosmological models for halos at redshift Zo — 0.05. Surprisingly, one 
observes that, contrary to what one might have expected, inhomogeneous models show 
less pronounced differences in both ps and A.v/2 than homogeneous models. Moreover, 
the inhomogeneous lines, in particular, do not reflect the anticipated Cyir trends at 
fixed mass. This occurs because ps and Ay/2 depend not only on Cvir (exactly on 
c^jr/(/n(cvii.) — Cvir/(1 + Cvir))) but also on A^i-^{zo) § which somewhat counteracts with 
Cvir- In effect, for inhomogeneous models, the non-linear overdensities at virialisation 
decrease as the concentrations increase. The homogeneous models, conversely, present 
differences among themselves that can be about a factor 1.5. In the next section 
we discuss the feasibility of the aforementioned techniques to actually discriminate 
cosmological models. 

Finally, the ability of a cluster to lens a background object, a galaxy or a quasar, 
depends critically on the halo profile and consequently on its concentration, therefore 
the influence of dark energy (and its homogeneity nature) on halo concentrations may 
be tested, albeit indirectly, through the observational statistics of multiply imaged 
quasars or giant arcs of galaxies. [12] have shown that the cross-section for quasar 
multiple imaging is increased for homogeneous dark energy models with w > —1, owing 
to halo concentrations being larger in those models. The expected increase on halo 
concentrations in models where dark energy is inhomogeneous may also produce an 
effect on the lensing cross-section. This remains to be investigated. 

5. Discussion and Conclusions 

We have studied halo properties in models with dynamical dark energy. This work 
extends upon previous studies in that we investigate halo properties as a function of a 
time-varying dark energy equation of state, covering four classes of potentials, and its 
homogeneity nature in the highly non-linear regime. 

The dynamics of the collapse of the halo is regulated by the form of the dark energy 
potential, the time evolution of its equation of state, the initial conditions of the field 
and its homogeneity properties in the highly non-linear regime. These have a direct 
influence on the turnaround, virialisation and collapse times, altering in consequence 
the (non-linear) density contrast and virial radius [15, 16]. 

As for the homogeneity nature of dark energy in the highly non-linear regime, we 
have examined two extreme scenarios: first, we assumed that the quintessence fleld 
does not exhibit fiuctuations on cluster scales and below - homogeneous fiuid; second, 

§ Ps has an additional dependence on Pcrit- 
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we supposed that the field inside the overdensity collapses along with the dark matter 
- inhomogeneous fluid. 

We have computed halo concentrations using the algorithm of [17] and have derived 
physical halo properties expected within this analytical treatment. We find that the 
[17] prescription displays, in general (the exception being the [29] model), higher 
halo concentrations in inhomogeneous dark energy models than in their homogeneous 
equivalents. The [36] inhomogeneous model is the model with the highest concentrations 
which are nearly a factor of two higher than the cosmological constant model or its 
homogeneous counterpart. For homogeneous dark energy models concentrations are 
larger as w increases owing to structures forming earlier. This is in agreement with 
previous findings from high-resolution numerical simulations which demonstrate that 
halos keep a memory in their central regions of the mean density of the universe at their 
characteristic formation epoch. 

In the two cases analysed (homogeneous and inhomogeneous models), halo 
concentrations seem to be controlled by both changes in formation epochs of the 
halo cores as well as by differing virialisation overdensities. While for inhomogeneous 
models, changes in the virialisation process constitute the most influential factor, for 
homogeneous models it is the differing formation histories that are most responsible for 
the changes in concentrations. 

Having determined the theoretical Cvir — Mvir relation, we then deduced the 
corresponding — and Ay/2 — Knax relations which represent physical measures 
that, unlike the Cvir — Afvir relation, may establish a more direct link with observations. 
Here, we note that the homogeneous models manifest more detectable differences in 
both ps and Ay/2 than the inhomogeneous models. Specially, the inhomogeneous curves 
do not register the scaling one might have expected based on the Cyir values at flxed 
mass. This arises because Ps and Ay/2 depend both on Cyir and on Ayij.{zo) which 
counterbalance each other. Nonetheless, it is intriguing why inhomogeneous models 
reveal visible changes in their collapse dynamics and do not show apparent differences 
in their central densities. This could be associated with the prescription used to calculate 
the characteristic collapse epoch which may not be applicable to inhomogeneous models. 
We notice that it is also important to realise that the free parameter = 28, whose 
meaning is in fact unknown, may not be appropriate. Clearly, these assumptions should 
not be undervalued, and require a test against N-body simulations. 

Studying the impact of dark energy on the density structure of dark matter halos 
surely represents an important step on our understanding of structure formation on 
those models. Unfortunately, halo properties do not seem to provide a (at present) 
robust way of probing dark energy. The Ay/2 — Knax and Ps — relations are plagued 
by observational scatter on the data that may (or may not) be associated with selection 
effects (cooling baryons in the former; hydrostatic equilibrium and isothermal hypothesis 
in the latter) and by the intrinsic, theoretical in nature, scatter in halo concentration 
values about the mean which is beheved to be due to a spread on the collapse histories 
of N-body simulated halos. Additionally, there is the implicit uncertainty on the the 
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value of (jg as a function of w which we fixed in this analysis (see discussion in [8]). 
Nevertheless, the effect of dark energy on the dark matter halo structure may be 
exploited indirectly through strong and weak lensing statistics. Another alternative 
path may be the time evolution of dark matter halo abundances [6]. 

While it is reassuring that the analytical recipe of [17] is successful in reproducing 
halo concentrations measured in numerical simulations performed in homogeneous dark 
energy models, more exphcitly for SUGRA and Ratra-Peebles potentials [7], it should 
be recognised that the same recipe has not been tested against simulations carried out in 
inhomogeneous models. Moreover, one should point out that the methods here used are 
based on the premise that the halo concentration-mass relation is governed by the ENS 
model. However, the ENS model is by no means the only viable and popular model. It 
is just one of a number of theoretical models for how the concentration - mass relation 
arises. See [18,19,44] for other examples of other possibilities. Furthermore, since the 
Dolag et al. paper [7] (the most recent paper advocating the ENS prescription) there 
have been a number of developments in studying the halo concentration - mass - redshift 
relation. Some of these studies produced results in agreement with ENS (e.g. Neto et 
al. [21]), others find a mass and redshift dependence that contradicts the ENS model 
(e.g. Maccio et al. [20] and Gao et al. [22]). At present the matter has not been settled, 
and it is not possible to identify one model that works better than all of the others. 
A Further investigation using and comparing other prescriptions will appear in a near 
future work [50]. 

Ideally, this work will motivate further studies using N-body simulations which, 
ultimately, will draw a firm conclusion on the interconnection between the dark matter 
halo structure and the underlying cosmology. 
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